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Abstract 

The Renormalization Group flow equations obtained by means of a proper time regulator are used 
to analyze the restoration of the discrete chiral symmetry at non-zero density and temperature in the 
Gross-Neveu model in d — 2 + 1 dimensions. The effects of the wave function renormalization of the 
auxiliary scalar field on the transition have been studied. The analysis is performed for a number of 
fermion flavors Nf = 12 and the limit of large Nf is also considered. The results are compared with 
those coming from lattice simulations. 
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In recent years there has been a renewed interest in the chiral symmetry restoration in QCD in 
connection with its possible experimental observation in the transition from the confined to the quark- 
gluon plasma phase. Lattice simulations, which are the most relevant tool to study such nonperturbative 
phenomena, when employed to investigate QCD at finite baryon density, must deal with the main problem 
of the appearance of complex terms in the Euclidean action due to the presence of the chemical potential 
and only recently new lattice techniques have been proposed P3G1, which overcome this problem, at least 
for not too large values of the chemical potential. On the other hand, the chiral symmetry breaking 
at finite temperature and density has been investigated in other models j3J |SJ d, which have a 
well defined action even for nonvanishing chemical potential and which, despite their simplicity, give a 
qualitative picture of the transition that could provide reasonable indications of this phenomenon in the 
full QCD framework. 

The simplest model is represented by the Gross-Neveu model jS], characterized by either discrete 
or continuous chiral symmetry, which in three spacetime dimensions shows a critical coupling g* that 
indicates the threshold for the occurrence of chiral symmetry breaking at zero temperature and density. 
This model has an interacting continuous limit with a spectrum that contains composite fermion and 
antifermion states analogous to baryons and mesons and its four fermion interaction can be regarded as 
an effective interaction for quarks at intermediate energies and at not very high temperatures where the 
effect of heavy mesons is still suppressed. Moreover in three dimensions the model is renormalizable order 
by order in the 1/N expansion around the critical coupling g* [jjj. 

In addition to lattice calculations, another nonperturbative tool, the Exact Renormalization Group 
(ERG) (see ^0] for a review of the subject and a more complete list of references) has been employed 
to study four fermion models, such as the Nambu-Jona-Lasinio model at finite temperature and density 
analyzed in [U], as well as the fixed point structure of the Gross-Neveu model in [12] . In this letter 
we make use of the Renormalization Group to investigate the specific problem of the phase diagram of 
the tridimensional Gross-Neveu model in order to explicitly test this approach through a comparison 
with the lattice analysis. Wc shall focus on a particular version of the Renormalization Group flow 
equation, based on the heat kernel representation of the effective potential, which has already been used 
in similar contexts, namely for analyzing the four-dimensional linear sigma model at finite temperature 
and densitv [T5IITllT5| . 

This latter version of the flow equation, often indicated as Proper Time Renormalization Group 
(PTRG) is typically obtained as an improvement of the one loop effective action by making use of the 
proper time representation \FS\ an d introducing a suitable multiplicative cutoff function as a regulator of 
the infrared modes [171 ITH| . The PTRG does not belong to the class of ERG flows which can be formally 
derived from the Green's functions generator without resorting to any truncation or approximation [191 
1201 121] . Nevertheless the PTRG flow equation has the property of preserving the symmetries of the 
theory and, in addition, it has a quite simple and manageable structure. Moreover the PTRG flow, 
although not able to fully reproduce perturbative expansions beyond the one loop order (221 7 provides 
excellent results when used to evaluate the critical properties such as the critical exponents of the three 
dimensional scalar theories at the nongaussian fixed point [281 124] . In particular the determination is 
optimized when one takes the sharp limit of the cutoff function on the proper time and it turns out to 
be much more accurate than the one corresponding to a smoother regulator. The PTRG also provides 
excellent determinations of the energy levels of the quantomechanical double well |2S]. These results are 
obtained at the first order in the derivative expansion of the effective action, i.e. by reducing the full flow 
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equation to two coupled differential equations for the potential and the wave function renormalization, 
and in our opinion this is a clear indication that, within this specific approximation, the PTRG flow is 
reliably accurate. According to this point of view it is interesting to test the PTRG flow in a different 
context which involves fermionic systems at finite temperature and density and in particular to compare 
the results coming from this approach with those obtained from the lattice investigations. 

Therefore in the following we focus on the Gross-Neveu model with Euclidean action 

= J d 3 x (iQ^ - 2^(^) 2 ) (1) 

where (j) = 7 M 9 P involves three gamma matrices (/i takes the three values: 1,2,4), rpj is a four component 
spinor and j is the flavor index : j = 1, . . .Nf. The action Q is invariant under the discrete chiral 
transformations ipj — > 7 tpj , ipj — > —i/j^j 5 . 

As it is well known 8 , instead of considering Eq. Q, it is helpful to introduce an auxiliary scalar 

field a and put the action in a bosonized form. In fact after adding to the term J d 3 x(l/2)(a + 
^2 



gipjipj/y/Nf) , which is quadratic in a and therefore does not modify the fermionic theory, one alterna- 
tively describes the model by means of the action 

S = Jd 3 x (^a 2 + ty-Wj + -J=a^^)j (2) 

The action (J2J) is to be taken as the bare action which fixes the ultraviolet boundary conditions of the 
renormalization group flow. In fact in the flow equations we consider a more general structure of the 
effective action which includes a scalar potential Uk(<r 2 ), the kinetic term for the scalar sector, properly 
normalized by the wave function renormalization Zk and also a fermionic wave function renormalization 



In Eq. Q we have explicitly introduced a dependence of the various terms on the running momentum 
scale k which parametrizes the PTRG flow. In fact the renormalization group equations describe their 
evolution as functions of k, starting from the ultraviolet boundary conditions given in Eq. J5J), down to 
the far infrared region, when all modes have been integrated out. So the ultraviolet boundary conditions, 
corresponding the large value k = A, are Za = 0, Z^ = 1, U\((T 2 ) = <j 2 /2, and a fixed value for g\. We 
shall limit ourselves to field independent wave function renormalizations Zk and and to a potential 
which is polynomial in a 2 . 

The shape of effective potential Uk(& 2 ) in the infrared region, i.e. at k ~ 0, does provide informa- 
tions about the chiral symmetry breaking and restoration, which are in fact related to the expectation 
value of the auxiliary scalar field a. If this expectation value is nonvanishing then the fermion becomes 
massive through the Yukawa coupling in © and the chiral symmetry of the bare action is spontaneously 
broken. As already mentioned, the mean field analysis indicates that at zero temperature and density 
the symmetry is broken only for couplings larger than a critical value, g > g* , whereas below g* the field 
a has zero expectation value and the fermion is massless, preserving chirality. When g > g* a sufficient 
increase of the temperature and (or) of the density is expected to restore the symmetry. 
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A saddle point expansion of Eq. yields the one loop effective action 

" S 2 S 
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The logarithms in Eq. can be expressed by means of the proper time representation 



ln(X) 



sX 



(4) 



(5) 



(0) 



and the infrared modes are properly regularized by means of a cutoff function fk, |18| . multiplicatively 
introduced into Eq. © by replacing e~ sX with e~ sX fk- Clearly this replacement introduces a depen- 
dence on the variable k in Eq. Then the PTRG equations are obtained by replacing everywhere in 
Eq. (@J the action with the k dependent ansatz Sk of Eq. Q and by taking the derivatives with respect 
to k of each side of Eq. Q . Typically fk is taken as a smooth heat kernel cutoff ^Hl 1131 1141 H3] , however 
as shown in |231 121) the PTRG equations are optimized when the sharp limit of fk is taken, or, more 
precisely when its derivative with respect to k ( which is the relevant quantity in the differential flow 
equations) is taken as a delta- function [201122 



kdkfk = ~2s S[s- 



k 2 



(7) 



Therefore in the following we take fk as implicitly defined in Eq. JJJ and the PTRG flow equation has a 
particularly simple form 



kd k S k = Tr exp 



Tr exp 



(Ztk) 2 



where we used a concise notation for the second derivatives of the action and Zk and Zf provide the 
correct normalization of the two exponential terms. The traces are extended to all the relevant degrees 
of freedom. 

The extension of the flow equation to finite temperature is done following the standard procedure |26j , 
i.e. by imposing periodic boundary conditions on the Euclidean time component with period fixed by the 
inverse temperature 1/T and consequently discretizing the corresponding component of the momentum. 
Therefore the integration on the component q± (the index 4 is taken in analogy to the one used for 
the gamma matrix 7 4 ) is to be converted into the infinite sum over the Matsubara frequencies, which 
are respectively 94 — > ui n = 2mtT for bosons and (74 — > v n = (2n + l)wT for fermions and n 6 Z, 
Jdq 4 ^2irTj: n . 

Finally, finite fermion density effects are accounted for by adding to the bare action the chemical 
potential fj, coupled to the density ip'ip. Then we can compute the purely fermionic contribution to Eq. 
||SJ in the momentum space 



2„2 



AiAk = (Z:) 2 ( q 2 + ^)+^ 2 +g 2 <T 



(9) 
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where Nf has been absorbed into the coupling constant g k y/Nf = g k . 

We are now able to extract from the full flow equation (JSJ the set of coupled equations for the various 
k dependent terms of the effective action. The equation for the potential U k is simply obtained by 
considering constant field configurations for a and vanishing fermion fields in Eq. (JHJ and, by performing 
the integral on the momentum, we get 

kd k u k = n b + n f (io) 

where the boson contribution, after performing the momentum integrals, is 



k 2 



M b 
Z k k 2 



(11) 



with M b = 2U k 
is 



Aa 2 U k (the dots indicate derivatives with respect to cr 2 ), and the fermion contribution 
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(12) 



Equation (|12l) shows the fermion contribution in a compact form. However the momentum integrals can 
be performed analytically which makes easier the numerical integration of the set of flow equations. 

In order to derive the equation for Z k one has to expand the exponentials in the right hand side of 
Eq. retaining the commutators of the coordinate and momentum dependent terms, and collecting 
the full coefficient of d^ad^a. This procedure is already discussed in and, specifically for the PTRG 
flow equation of Z k in scalar theories, in [2311241 122] . Moreover here we are only interested in the leading 
corrections to the flow of the potential U k and therefore, instead of solving the problem by taking Z k , as 
well as Zf and g k , as arbitrary functions of the fields, we we limit ourselves to consider the particular 
value of these parameters when computed at the minimum of the potential at each value of k, a = a k 
and for vanishing fermion fields. The lowest order equation for Z k is therefore: 



kd k Z k = - 

TN f g 2 
(Ztk) 2 n 
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(13) 



The evolution equation of Zf is obtained by taking the coefficient of ipjftipj after expanding the 



exponentials in the right hand side of Eq. (|H| 
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(14) 
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and the one for g k by selecting the coefficient of aipjipji 

where Mt and Mj, in Eqs. ljT3|) . itTH . ((T3|) are evaluated at a = Wk and W is defined as 

W = zf( q 2 + v 2 ) + g 2 W 2 ~ f , 2 (16) 

However, in the following we shall not solve the full potential equation Eq. (|10l) which is a partial 
differential equation, but rather we treat a set of ordinary differential equations for the various coefficients 
of a truncated polynomial expansion of the potential. In fact, according to |14j such truncations converge 
quite rapidly and we choose to retain powers of the field up to er 8 in the expansion. Therefore we 
parametrize the potential, expanded around the local minimum Wk , in the following way 

I 4 
U k (<j 2 ) = w , fc + -W2,fc(o- 2 - CTfe 2 ) + ^u 2 i : k(cr 2 ~ek 2 ) % (17) 

i=2 

The parametrization in Eq. p7|l contains five k dependent parameters. In fact when the minimum of the 
potential is located at Wk — 0, then the independent parameters are uo./c, W2,fe, ^4,fcj u 6,k, us,k- When 
instead the minimum is at a nonvanishing Wk, one has Uk(W 2 ) = and then U2,fc = identically and the k 
dependent parameters are still five with v,2,k replaced by Wk- For our potential at k = A we have 1*2. a = 1 
and the minimum at Wa — and therefore we start our flow by explicitly considering the evolution of 1*2, 
Then, when k is lowered, U2,k diminishes and (only for suitable values of T and /1) eventually vanishes 
at k ( with < k < A). So, at k = k, we replace the flow equation for 112, k with the one for Wk, with the 
initial condition W-r = and require continuity for the other parameters tio,fe, U4,fc, u&,k, Ut8,k- 

The flow equations for the parameters in Eq. I|17|l are directly obtained from Eq. I|1U|) by repeated 
differentiations with respect to a 2 (we shall not display them explicitly) and the one for W k is instead 
derived from the minimum condition, 

(2W k Uk(W 2 k ))kd k Wk = -kd k U k (W 2 k ) (18) 

Incidentally we note that ito.fc never appears in the other equations and therefore this parameter does 
not affect the evolution of the others. 

As a first step in the analysis of the flow equations, we briefly review the well known case Nf — > 00, 
that can be treated analytically. In fact, since the minimum of the potential increases as yNf while 
the generated mass of the fermion does not depend on Nt (we recall that in our definitions the coupling 
gk is proportional to 1/WNf ) it is convenient to rescale the fields as a — > y/Nfcr, iji — * y/Nfip as well 
as the potential Uk — > NfUk so that the action, expressed in terms of these new variables and of gk, is 
homogeneously rescaled by the factor Nf. It is easy to realize that after these changes the right hand 
sides of the flow equations for and gk are suppressed in the limit Nf — > 00 and these parameters 
remain fixed to their initial value at k — A and in particular Zf = 1. Also, in Eq. I jlOp the term 
is suppressed and the equation for the potential depends only on the fermionic contribution which does 
not contain any dependence on the potential itself or on its derivatives or on the parameter Z k - As a 
consequence, the equations for Zk and for Uk decouple. We neglect the former and focus on the latter, 



Z k (q z +^ n ) + M b 



Z k k 2 



W 



(15) 
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which is relevant for the symmetry breaking problem, and we note that in this limit Eq. (|10fl is no longer 
a differential equation and can be simply integrated in k from zero to infinity with the help of Eq. © 
and by recalling the boundary condition for the potential in the ultraviolet region assigned in Eq. J2J. 
By extending the k integration up to infinity, we have implicitly removed the ultraviolet scale A and this 
will introduce some divergences in our computation which at some point need to be regulated. However 
in this way we can more easily make contact with the standard calculations of the large Nf case (see e.g. 
|28|V In order to determine the extrema of the potential it is convenient to focus on its derivative with 
respect to a and, after performing the sum over the Matsubara frequencies PS], we get ( we neglect for 
simplicity the subscript k for the coupling g since here, as noticed above, it is k independent ) 



dUk=G 
da 



9 2 « 



(2tt)2 v /^T 



gfo 2 



tanh 
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2T 



(19) 



In the particular case fj, — T — 0, if we regulate in Eq. I|19f) the divergent integral on the momentum q 
with the ultraviolet cutoff A, we get the gap equation that defines the fermion mass at zero temperature 
and chemical potential (28] 

i= A-ga = A- TO/ (20) 
9 

If we consider T = and finite chemical potential we see that the hyperbolic tangents in Eq. I|19fl in 
the limit of vanishing temperature become step functions and, by making use of the gap equation (|2U|) 
to get rid of the regulator A, it is easy to check that taking /i < to/ or \i > rrif yields two different 
behaviors. In fact in the former case the same gap equation as for /x = is obtained, in the latter case 
we get the extremum condition g 2 c(/i — rrif) — that has the unique solution a = 0. Therefore at T = 
a first order transition is observed when the chemical potential is equal to to/. 

Finally when also T^O, again with the help of Eq. l]2Up. we get the gap equation 



cosh 



T 



,(m s /T) 



cosh- 



T 







(21) 



For nonvanishing T, Eq. 1)21(1 . in addition to a = 0, has another solution only for sufficiently small 
chemical potential \i < Ji(T) ( with ~p(T) ~ m/ for small T ) . This solution tends to zero for /i — > Jl(T) 
from below and disappears for /i > ~p(T) where the only extremum is a — 0. Then the transition in this 
case is second order. Only at T = there is a singular behavior that leads to a first order transition [5]. 

Let us now consider the finite Nf case and, in order to test the PTRG we shall compare our analysis 
with the lattice results obtained in [S] where Nf = 12 is taken to observe the leading corrections to the 
Nf — * oo case. Therefore in the following we select this particular value of Nf and solve numerically 
our flow equations. The flow starts at the scale k = A and ends at k = and we varied A in the range 
10 3 — 10 5 to test the insensitivity of the final output of the dimcnsionless parameters to the specific value 
of the ultraviolet scale. The initial values of the parameters at k — A are those indicated below Eq. © 
and it must be noted that it2,A = 1, independently of the particular value of A and this implies that 
U2,k must be dimensionless. Consequently all the canonical dimensions of the other parameters must be 
arranged according to this point. 

The most simple version of the Sow equations is obtained by a further reduction of the number of 
running parameters and, specifically, by neglecting the flow of the renormalization of the scalar field. 
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Since in our problem the bare value of this parameter is = 0, this approximation corresponds to 
have Zk = at any scale k and therefore to totally neglect the kinetic contribution of the scalar field. 
Moreover, as it can be seen from Eqs. (|14|l and (|15|l . in this case also Z% and gk become ^-independent 
and stay fixed to their initial value and in Eq. 1)10(1 . fib = 0. When Zk is turned on and the full set of flow 
equations is considered, still Z^ and gk do not show any significant change from their initial values. On 
the contrary Zk has a not negligible flow and when the scale k goes from the ultraviolet to the infrared 
region, this parameter grows from zero to a finite value, and this was already observed in models of 
fermions coupled to scalars [111 I29| . In general the effects of the scalar renormalization are significant, as 
it happens for instance in the quantitative determination of the critical exponents of the scalar theory or 
of some energy levels in quantum mechanics |24l I25j and, also in this problem they will turn out to be 
relevant. 

At fi = T = and with Nf — 12 (which is the value of Nf adopted in we find a critical 
value of the coupling that separates the symmetric and the broken phase, which, in cutoff units, is 
l/(g\VA) 2 = 0.178(2). In |S] the inverse squared critical coupling, expressed in different units, is 0.98 
and, by addressing the difference between these two determinations to the two different scales employed 
to normalize the dimensionful coupling, it follows that the ratio of these two numbers directly yields the 
ratio of the two scales. Then, since the value of the inverse square coupling taken in [5] to study the 
restoration of the chiral symmetry at finite T and fi is 0.75, we can easily convert it in units of our cutoff 
A to get the input value of the coupling in our flow equations. In our scale it is l/(gAV / A) 2 = 0.136. At 
this value of the coupling and again at fi = T = we have determined the fermion mass which, according 
to the criterion chosen in [S], will be the parameter used to normalize the dimensionful quantities in our 
analysis at finite fi and T. In particular we found rrif /A = 0.146 and, when Zk = 0, rrif /A = 0.145. 

The phase diagram found by means of the flow equations is shown in Fig. 1. The solid line shows the 
transition at Nf — ► oo as obtained from Eq. I|21|) . Except for the points at fi — where the temperature 
is raised in order to observe the symmetry restoration, all the other points are obtained by keeping 
a fixed temperature and changing the chemical potential fi to induce the transition and the values of 
the temperature considered are: T/rrif = 0, 0.24, 0.29, 0.36, 0.48. For the circles at fi = and at 
T/mf = 0.48 the errors are small and they are not displayed in the figure. Notice that the circle at 
/i = is almost coincident with the lattice determination of [Jj which corresponds to the star (together 
with its error bar) on the T axis in Fig. 1. As expected the circles are in very good agreement with the 
mean field line and it is necessary to go beyond the local potential approximation to improve the results. 
The inclusion of the flow of Zk induces a significant change, reducing the values of T and fi at which 
the transition occurs. This effect, qualitatively correct, is somehow large and the points turn out to be 
systematically lower than the lattice results. 

Also, some uncertainties affect our determinations of the transition points in the fi — T plane as it is 
shown by the error bars associated to the points in Fig. 1. In fact, for values of fi below these bars ( or, 
in the case of the transition at fi — 0, for values of T below the error bar) we find the absolute minimum 
of the potential at a nonvanishing vacuum expectation of the field whereas, for values above these bars 
the minimum corresponds to a vanishing vacuum expectation of the field and chiral symmetry is thus 
restored. For values of fi (and T) within the bars it is difficult to determine the nature of the transition, 
especially if the flow of Zk is included and therefore, by taking the most conservative point of view, we 
explicitly indicated in Fig.l the ranges of fi and T where these uncertainties are present. 
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As an example, in Fig. 2 ~ak=ol m f is displayed as a function of T/mj. The solid line is obtained for 
Zk = and it shows a continuous transition and the dashed line corresponds to Zk ^ 0. The latter curve 
has a discontinuous fall to Uk = which should indicate a first order transition. But in the temperature 
range indicated by the error bar on the £-axis of Fig. 2 (which is the same error bar shown for the 
corresponding point in Fig. 1) we observe in our flow equations the appearance, at small finite k, of a 
non- vanishing minimum of the potential which however, still at finite k, vanishes. Then, for temperatures 
above the error bar in Fig. 2 the minimum of the potential corresponds to at — at any value of k. 
Because of this effect that , incidentally, has already been observed in J3| > we believe that it is not safe 
to conclude that the inclusion of the flow of Zk has changed the order of the transition which moreover 
would be in contrast with the expectation of a second order transition. Rather, we would argue that it 
is more likely that a different truncation in our flow equations, which includes more operators, like field 
dependent contributions to the scalar wave function renormalization, could restore the smooth behavior 
obtained when Zk = and the dashed line in Fig. 2 would continuously reach zero as it happens for the 
other one. 

In conclusion the PTRG flow equations provide a good qualitative picture of the transition at finite 
temperature and chemical potential, and in particular they have the correct behavior in the limit Nf — > oo. 
Also Fig. 1 shows that the role of Zk is certainly relevant in the determination of the critical line and a 
reasonable comparison with the phase diagram obtained from the lattice simulations in J5]j is possible. 
As a further check we have computed, following [5], the exponent v for the critical behavior around the 
coupling g* at T = 0, and we found v — 1.06(4) when Z% = and v = 1.19(30) for the full set of 
equations, which has to be compared to the value v = 1.05(10) found in 5 . Again the inclusion of Zk 
gives a correction larger that the one expected. On the other hand, the analysis of the order of the phase 
transition at this level of approximation is not satisfactory, although some expected features are recovered. 
In our opinion the quantitative discrepancies as well as the various uncertainties discussed have to be 
addressed to the specific truncations and approximations made on the flow equations. We expect that 
the inclusion of field dependent terms in the scalar wave function renormalization (as it has been already 
pointed out in |15| in a different context ) could significantly improve the analysis and therefore this turns 
out to be an essential ingredient to carry out an accurate analysis of phenomenologically realistic models. 
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Figure 1 
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Figure 1: Critical line on the plane T — fi. The solid line corresponds to the transition for Nj — > oo. 
The circles are obtained from the flow equations with Zk = and the diamonds for Zk =/= 0. The stars 
correspond the lattice results as quoted in T and /j, are expressed in units of the fermion mass m/. 
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